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Abstract 

We report a theoretical scheme using a B-sphne basis set to improve the poor computational 
accuracy of circular Rydberg states of hydrogen atoms in the intermediate magnetic field. This 
scheme can produce high accuracy energy levels and valid for an arbitrary magnetic field. Energy 
levels of hydrogen are presented for circular Rydberg states with azimuthal quantum numbers \m\ 
= 10 - 70 as a function of magnetic field strengths ranging from zero to 2.35 x 10^ T. The variation 
of spatial distributions of electron probability densities with magnetic field strengths is discussed 
and competition between Coulomb and magnetic interactions is illustrated. 

PACS numbers: 32.60.+i, 03.65. Sq 
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I. INTRODUCTION 



Circular states of Rydberg atoms are states with \m\ = n - 1, where m and n are the 
azimuthal and principal quantum numbers, respectively. The word "circular" here is used 
because the classical orbits corresponding to these quantum states are circular. Such states 
have drawn considerable interests since the early development stage of quantum theory. In 
modern physics, circular Rydberg states play a crucial role in understanding the correspon- 
dence limit between classical and quantum physics. 

Rydberg atoms in circular states possess some special properties, such as the largest mag- 



netic moment, the smaf 
collision cross sections 



est Stark effect, the longest radiative lifetime, and highly anisotropic 



l2|. 



mains, including cavity QED 



hese properties make them applicable in many research do- 
3j, quantum computation and information processing 



measurements of fundamental physics constants ultrafast dynamics of Rydberg atoms 
Q], and so on. The earliest experiment for circular state populations was made by Hulet 
and Kleppner In that experiment, the circular states of lithium atoms with |m| = 17, 18 
were prepared using an adiabatic microwave transfer technique. Soon afterwards, the same 
technique was applied to the other two alkali metal atoms, cesium [8| and rubidium [9|. 
The external electric and magnetic fields have been adopted to control adiabatic microwave 
transfer and to split ±m degeneracy, respectively {sl, A theoretical scheme using crossed 
electric and magnetic fields to populate circular states was proposed by Delende and Gay 



ll| and realized experimentally by Hare et al. [12] . The creation of circular states of hydro- 



gen atoms was based on both adiabatic microwave transfer and the crossed fields method 
2j. Very recently, Maeda et al. 13| successfully made long-lived, nondispersing Bohr wave 
packets starting from Li Rydberg atoms in circular states. Obviously interests in circular 
states have been enhanced in recent years. 

Besides fundemental physics interests, investigation of magnetized H atoms is of applica- 
tion importance in astronomy and astrophysics. To understand the behavior of magnetized 
white dwarf stars (10^ — 10^ T) and neutron stars (10'' — 10^ T), the accurate knowledge 
of the physical properties of atoms in strong magnetic fields is essential. 

Much attention has been paid to experiment and astronomical observation, but corre- 
sponding computations of circular Rydberg states of atoms in external fields are relatively 
rare. Wunner et al. [l^ solved the Schodinger equations for circular states of magnetized 
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hydrogen atoms. For low fields, they expanded wave functions in terms of a complete, or- 
thonormal set of functions in the radial direction and spherical harmonics in the angular 
direction, while for high fields (< 10'' T), they employed the adiabatic approximation. In the 
intermediate field range, they smoothly jointed the two-side results. Seemingly, the method 
of Wunner et al. may effectively handle any magnetic field from low to high. However, that 
is not true. It is widely known that the adiabatic approximation is valid only for extremely 
high magnetic fields, and interpolation in the intermediate field range may produce poor 
results because of the used energies stemming from the adiabatic approximation. Here it is 
necessary to point out that intermediate magnetic field used in literature is a very ambiguous 
concept. Strictly speaking in physics, the "intermediate magnetic field" means a magnetic 
field when Coulomb and magnetic interactions are of a comparable magnitude. In other 
words, also it is closely relevant to the energy levels besides magnetic fields. 

Liu et al. 15(| performed calculations of circular Rydberg states of H atoms in magnetic 



field strengths ranging from zero to 70.5 T, based on a B-spline expansion method, in 
which wave functions are expanded in terms of B-spline functions in both the radial and 
angular directions. However, their method is limited to low magnetic fields. It should 
be emphasized, again, that the ambiguous concept of the so-called low magnetic fields is 
)eing used here. The disadvantages of their method have been analyzed recently in Ref. 



16|. The accuracy of their energy levels decreases with magnetic fields B when B > 47 



T. Germann et al. [17i] developed a dimensional perturbation theory to evaluate circular 
Rydberg states of magnetized atomic hydrogen with an extreme high accuracy. The theory 
is applicable to entire range of magnetic field strengths, but is limited to high azimuthal 
quantum number (m ^ 1). Later, Watson and coworker 18|] extended their method to 
include low-m situations. However, not in all field strengths, circular states may be treated 
because of the limitation of this theory. 

This paper focuses on the problem of circular Rydberg states in the intermediate mag- 
netic fields. First of all, we develop a method to evaluate circular Rydberg states of atomic 
hydrogen in a strong magnetic field. In cylindrical coordinates, the wave functions are ex- 
panded in terms of Landau states in the p direction and the B-spline basis in the z direction. 
This method is combined together with a recently reported finite-basis-set technique 16| . 
which is valid for weak fields. We will show that the combination scheme is a high accuracy 
tool to compute of circular Rydberg states of atomic hydrogen in the intermediate magnetic 
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fields, and valid for an arbitrary magnetic field. 



II. THEORETICAL METHOD 

The nonrelativistic Hamiltonian of hydrogen atoms on a uniform magnetic field along the 
z axis reads in atomic units 

where 7 is the magnetic field strength in units of Bq 2.35 x 10^ T, and are the z 
components of the orbital and spin angular momenta, respectively, the third term linear 
in 7 is the paramagnetic potential, and the fourth term quadratic in 7 is the diamagnetic 
potential. It is easily seen that H commutes with (corresponding quantum number m) as 
well as Sz- Depending on relative magnitudes of Coulomb and diamagnetic potentials, the 
symmetry of the system changes. The system is more spherically symmetric if the Coulomb 
potential is dominant, and more cylindrically symmetric if the diamagnetic potential is 
dominant. 

For strong magnetic fields, we solve the Schrodinger equation in cylindrical coordinates 
by expanding wave functions in terms of Landau states in the p direction and the B-spline 
basis in the z direction. The wave function is expressed as 

^{p,z,<f)) = J2C^nB,{z)n^{p,<f)), (2) 

i,n 

where TZn{p,4>) is the normalized wave function for the Landau state with n > 0, given by 

H 



,(l™l+l)/2 



7^„(p,0) = ^„(p)-= = ^^^^^ /'!!iL^ f e-'^'^/'p\-KF,i-n, \m\ + 1, pV2)' 



271 m\ V 2lHn! r . , , , ,r ,/ 

(3) 

with being the confluent hypergeometric function, and Bi{z) are the B-spline functions. 



The B-spline functions along with their properties have been out 
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ined in Ref. [l^, and 

Therefore, here it 



a complete description can be also found in the book of de Boor 
is unnecessary to give any elaboration once again. 7^„(p, 0) in Eq. ([3]) should have been 
denoted as 7^5^(p,0), but the superscript m is dropped for simplification. Considering m is 
a good quantum number, doing so should not give rise to any confusion. Substituting Eq. 



([2]) into the Schrodinger equation with the Hamiltonian ([T]), and projecting onto the basis 
BiTZn yield the following matrix equation 

HC = EAfC, (4) 

where the Hamiltonian matrix elements are of the form 

Hi'n'M = -2^n'n J Bi' (z) —Bi{z)dz + U-fSn'n J Bi> (z) Bi{z)dz 

+ / B,{z)Vn'n{z)Bi{z)dz, (5) 

Jo 

and TV is the overlap matrix with elements 

P ^max 

M'n';in = ^n'n Bi, (z) Bi{z)dz. (6) 

Jo 

The upper limits of the integration, Zmax, in Eqs. fl5ll6l) denotes the maximum of coordinate 
z. It changes with energy levels and magnetic field strengths. Zmax should be taken to 
be large enough to obtain stable eigenvalues in calculations. The overlap matrix A/" is a 
positive definite real symmetric matrix. Its existence is attributed to the B-spline basis 
nonorthogonality. The matrix element Vn'n{z) in Eq. ([5]) is given by 2l| 



Vn'n{z) = - I Mn'ip)^==^nip)pdp 







2 -f^-^ K!^!K + H)!(n+H)!]V^ 

2^ (n' - nV(n - ^Vdrr^l ^ nVn\ -^^"^ ^J, + |m| , , ^^J 



with 



I{a, l3; z) = J (sina;)^'^(cosx)^" exp 



7z^cos^ X 



dx. (8) 



f7r/2 

Gin ^^^^ rrtriG 'T^^^^ fivn 

2sin^ X 

The adaptive Gauss-Kronrod quadrature and the Gaussian quadrature are, respectively, 
used to evaluate the integral of Eq. ([8]) numerically, and all the matrix elements of H and 
M . For states of even and odd parity, we perform calculations at the whole space of z with 
—Zmax < -2 < Zmax and the half space of z with < z < Zmax-, respectively. The combination 
of such a choice and properties of B-spline functions make it easy to enforce boundary 
conditions. That is to simply remove i = 1,N terms from the summation over i in Eq. ([2D, 
namely -81(2;) and B]\f{z), where N represents the number of B-spline functions. Once the 
calculations of the matrix elements are finished, standard routines for matrix diagonalization 
can be used to obtain the eigenvalues and eigenvectors numerically. 



It should be pointed out that this method is vahd only for strong fields, because the 
basis size required to obtain stable convergence becomes larger and large with decreasing 
ield strengths. In the situation of weak fields, we will adopt the finite-basis-set technique 



16|, in which the Schrodinger equation is solved in spherical coordinates by expanding wave 
functions in terms of a B-spline basis in the radial direction and spherical harmonics in 
the angular direction. The finite-basis-set technique was developed recently in order to 
calculate ground and lowly excited states for magnetized hydrogen atoms. Its details have 
been presented and therefore are omitted here. We will combine the current method and the 
finite-basis-set technique for computations of circular Rydberg states of atomic hydrogen for 
arbitrary magnetic field strengths. 

III. RESULTS AND DISCUSSION 

Circular Rydberg states with |m| = 10 — 70 were calculated in a wide field strength 
range. Only the eigenvalue problem for states with negative m is solved, since eigenvalues 
for states with positive m are easily obtained by an energy shift from the corresponding 



states of negative m 14| . The maximum number and order of B-spline functions were taken 
to be = 100 and k = 11, respectively. A crucial procedure in the current calculation is 
to suitably select the knot sequence of B-spline functions and maxima of space coordinates, 
Zmax and Tmax, whcrc Tmax represents the maximum of radial coordinates r. For calculations 
in cylindrical coordinates with the current method, z = {0, Zmax} was divided into the inner 
and outer region, and in each region an equidistant mesh was taken. Furthermore, the mesh 
size in the inner region was selected to be equal to twice the mesh size in the outer region. 
Such a choice is expected to characterize the behavior of circular electron motion in the 
combined Coulomb and magnetic fields quite well. For calculations in spherical coordinates 
with the finite-basis-set technique, an exponential mesh was adopted, following Ref. [l^ 
(see this reference for details). Energy levels were optimized by adjusting the distribution 
of the knots, and z^ax (or r^a^.). 

For any given state, there does exist a magnetic field range with strong competition 
between Coulomb and magnetic interactions. In this range, the solution of the Schrodinger 
equation is very difficult. Either a gap of fields, where energy levels can not be directly 
calculated because of the limitation of the used method, remain left or the gap is smoothly 
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jointed to the results from the two sides of the gap (see, i.e., Ref. [ij] and references therein). 
Our method, however, is very effective in such a field range. Table I lists eigenvalues for 
circular states with m = — 10 and — 70 in a magnetic field with strengths ranging from 
0.0001 to 0.1 a.u., calculated in the current method and the finite-basis-set technique. The 
uncertainty of each eigenvalue is ± 1 in the last digit, except for those in a range from 0.0001 
to 0.002 a.u. in the second column and from 0.02 to 0.1 a.u. in the third column. The two 
methods are in excellent agreement. 

The numbers of spherical harmonics and Landau states, which are used to expand wave 
functions, are also given in this table. The numbers reflect symmetry of the system at a 
given state and field strength because of the competition between Coulomb and magnetic 
interactions. For example, let us look at the circular state with m = — 10. At 7 = 0.0001 
a.u., the system is more spherically symmetric (the number of spherical harmonics is 2), as 
the Coulomb interaction is dominant, while the magnetic interaction is relatively small, but 
at 7 = 0.1 a.u., the symmetry is completely broken (the number of spherical harmonics and 
Landau states is close), as the two kinds of interactions become comparable. Now, let us 
turn to the m = — 70 state. At 7 = 0.1 a.u., the system is more cylindrically symmetric (the 
number of Landau states is 5), as the magnetic interaction is dominant, while the Coulomb 
interaction is relatively small, but at 7 = 0.0001 a.u., the symmetry is completely broken, 
because the two interactions become comparable. 

It is time saving to adopt the different method for systems with more spherical and more 
cylindrical symmetries. Totally, the method developed in this paper is not suitable for weak 
magnetic fields, while the finite-basis-set technique is not suitable for strong fields. However, 
both methods remain valid for intermediate fields, and should be in excellent agreement. 
Eigenvalues are calculated and listed in Table II for circular Rydberg states with m from — 
20 to — 60 in a magnetic field range from zero to 10^ a.u. We first used the finite-basis-set 
technique when the magnetic field is small and then switched to the current method when 
it becomes larger than some value. 

We made a comparison between thepresent eigenvalues for the circular state with m = 



— 24 with available theoretical data 



14 



15 



171], and excellent agreement is found. The 



compared data were not listed here as doing so is obviously superfluous. Restricted by the 
used spherical coordinate system, the method of Liu et al. [l5| is ineffective for higher fields. 
One may see that their accuracy decreases with increasing field strengths. At 7 = 0.0003 
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a.u., our energy agrees very well with the highly accurate result of Germann et al. [iTJ. 
Our calculation gives —5.315 899 713 035 a.u. As previously mentioned, here we point out 
again that not in all field strengths circular states may be calculated with the dimensional 
)erturbation theory of Germann et al. because of the limitation of this theory (also see Ref. 



18|). 

Gontour plots of probability density distributions of the electron in a circular state with 
m = — 24 are drawn in Fig. 1 at several field strengths. This figure illustrates the change of 
probability density distributions of the circular Rydberg electron with magnetic fields. It is 
seen that the electron cloud becomes squeezed toward the 2;-axis with increasing magnetic 
fields. At the field-free case, the probability density widely distributes near p = 620 a.u., but 
the electron cloud shrinks in a very narrow region near p = 1 a.u. when the field strength 
increases up to 7 = 100 a.u. Furthermore, it is shown that the peak value of probability 
densities at 7 = 100 a.u. is 5 orders of magnitude larger than that at 7 = a.u. (see the 
color bars of the figure). 

To illustrate competition between Goulomb and magnetic interactions, probability density 
distributions V{p, z) of the electron in a non-circular state n^m = 7q_4 are drawn in Fig. 2 
at several magnetic field strengths. The readers are referred to Ref. 16| for the definition 
of state nirn- Also we calculated V{p, z) using the analytical formula of wave functions for 
the field-free case and found excellent agreement with the current numerical results. At 7 = 
0.003 a.u., the figure shows the system is more spherically symmetric. Only six spherical 
partial waves are needed to obtain the energy level convergent up to 10 digits of precision. 
From 0.006 a.u., the outermost and intermediate electron clouds are separated because the 
magnetic interaction is becoming significant. At 7 = 0.1, the Lorenz force is dominant. 
To reach the same convergence as 16 Landau states, we contained the 38 spherical partial 
waves. 



IV. SUMMARY AND CONCLUSION 

A theoretical scheme using a B-spline basis set has been proposed to compute circu- 
lar Rydberg states of hydrogen atoms in the intermediate magnetic field. We solved the 
Schrodinger equation in the spherical and cylindrical coordinate, corresponding to weak and 
strong magnetic fields, respectively. The scheme is valid for an arbitrary magnetic field, and 
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can produce high accuracy energy levels (Note: the high accuracy here indicates that in 
the theoretical framework defined). We showed how to save computational time as well as 
to obtain a high accuracy. The energy levels of atomic hydrogen are presented for circular 

Rydbcrg states with azimuthal quantum numbers \m\ = 10 - 70 as a function of magnetic 
field strengths ranging from zero to 2.35 x 10^ T. The competition between Coulomb and 
magnetic interactions is discussed and illustrated. 
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TABLE I: Eigenvalues for circular Rydberg states with m = — 10 and — 70 in units of Hartree in an intermediate magnetic field strength, 
calculated with the finite-basis-set technique (Sph.) and the method developed in this paper (Lan.). 



m = - 10 



m 
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B] = Ax 10^ 

The number in parenthesis gives the number of spherical harmonics. 
The number in parenthesis gives the number of Landau states. 



TABLE II: Eigenvalues for circular Rydberg states with m from —20 to —60 in units of Hartree and as a function of magnetic field strengths 
ranging from 7 = to 10^. The uncertainty of each eigenvalue is it 1 in the last digit. 
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Figure Captions 



Fig.l (Color online) Contour plots of probability densities of the electron in a circular state 
with m = —24, integrated over the angular variable of cylindrical coordinates (p, 0). 
The magnetic field is along the z-axis. The total wave function is normalized such 
that / z, (j))\^pdpdzd(t)=l. Note: The data aspect ratio is taken to be 1:1 between 
z and p. 

Fig.2 (Color online) Probability densities of the electron in a non-circular state Iq-a, inte- 
grated over the angular variable of cyhndrical coordinates {p,z,(j)). The direction of 
magnetic fields and the normalization of wave function are the same as those in Fig.l. 
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